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Abstract 

We outline a relationship between conformal field theories and spectral prob- 
lems of ordinary differential equations, and discuss its generalisation to mod- 
els related to classical Lie algebras. 



1 Introduction 



The ODE/IM correspondence [U EJ [3lll] has estabhshed a hnk between two di- 
mensional conformal field theory (CFT) and generalised spectral problems in or- 
dinary differential and pseudo-differential equations. It is based on an equivalence 
between transfer matrix eigenvalues O |6] and Baxter Q-functions in integrable 
models (IMs), and spectral determinants [Tj |8] of ordinary differential equations 
(ODEs). 

In statistical mechanics, the transfer matrix and its largest eigenvalue - de- 
noted by T in the following - are central objects. For example, consider the 
six-vertex model defined on a square lattice with N columns and N' rows; T can 
be written in terms of an auxiliary entire function Q through the so-called Bax- 
ter TQ relation. Up to an overall constant, Q is completely determined by the 
knowledge of the positions of its zeros, the Bethe roots, which are constrained by 
the Bethe ansatz equations (BAE). Subject to some qualitative information on 
the positions of the Bethe roots, easily deduced by studying systems with small 
size, the Bethe ansatz leads to a unique set of ground-state roots. In the A^' — > oo 
limit the free energy per site / is simply related to T by 

/~-— InT. (1.1) 

In [51 [6], Bazhanov, Lukyanov and Zamolodchikov showed how to adapt the 
same techniques directly to the conformal field theory (CFT) limit of the six- 
vertex model. In this setting, we consider the conformal field theory with Virasoro 
central charge c = 1 corresponding to the continuum limit of the six- vertex model, 
defined on an infinitely-long strip with twisted boundary conditions along the 
finite size direction. The largest transfer matrix eigenvalue T depends on three 
independent parameters: the (rescaled) spectral parameter v, the anisotropy rj 
and the twist (j). Defining E, M, I, uj, 0, through the following relations 

9„ M ^ ^ 9«,r , (21 + 1)tt 

E = e^'', 7]=-— , uj = e^, Jl = w2M cP= \; ' 1.2) 

''2M-hl ' ^ 2M + 2 ^ ' 

the resulting TQ relation is 

21+1 21+1 

T{E, I, M)Q{E, I, M) = UJ ~Q{nE, I, M) + uj~Q{Q.~^E, I, M) . (1.3) 

The Baxter function Q for this largest eigenvalue is fixed by demanding entirety 
of both T and Q, and reality, positivity and 'extreme packing' for I > —1/2 of 
the set {Ei} of zeros of Q. The BAE follow from the entirety of T and Q via 

Q{E,) = ^ nEME,) = ^ = -^^'-^^ • (1.4) 
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Surprisingly, equations (jl.3p and (jl.4p also emerge from an apparently unrelated 
context: the study of particular spectral problems for the following differential 
equation 

with X and E possibly complex. To see the emergence of ()1.4p from (jl.Sp . we 
start from the unique solution 'i/;(x, -B, /) of (jl.5p on the punctured complex plane 
X G C \ {0} which has the asymptotic 

V^~x-*^/2exp(-^x^^+i) , (M>1) (1.6) 
as |x| — > CX3 in any closed sector contained in the sector | argx | < 2K1+2 ' 

This 

solution is entire in E and x . From ijj we introduce a family of solutions to (jl.5p 
using the 'Sibuya trick' (also known as 'Symanzik rescaling'): 

V'fc = 'ilj{uj''x, n^E, I) . (1.7) 

In (jl.7p . A; takes integer values; any pair {ipkjipk+i} constitutes a basis of solutions 
to (jl.Sp . An alternative way to characterize a solution to (jl.Sp is through its 
behaviour near the origin x = 0. The indicial equation is 

(A-1-/)(A + /) = , (1.8) 

and correspondingly we can define two (generally) independent solutions 

X+(x, E) = x{x, E, I) ~ x'+i + 0(x'+3) , (1.9) 

and x~ixj -^) = xi^j — / — 1), which transform trivially under Symanzik rescal- 
ing as 

X+ = X+(^'a;,J]^ii;)=a;('+i)V(x,i?) . (1.10) 

The trick is now to rewrite xt = X^(^)-^) respectively in terms of the basis 
{V'cV'i} and {iI^-i.iIjq}: 

2ixt = uJ-'-'^Q{nE)i:o-Q{E)Lo-h^^ (1.11) 
2ixt = 2iJ+W = u2Q{E)i,_^-J+hQ{n-^E)iljo (1.12) 

where the coefficients has been fixed by consistency among (jl.lip . ()1.12p and 
(fLTOll and 

Q{E,l) = W[i,o,xt\- (1-13) 

Here <?] = — denotes the Wronskian of / and g. Taking the ratio 
(jl.lip / (jl.l2p evaluated at a zero E=Ei of Q leads immediately to the Bethe 
ansatz equations ()1.4p without the need to introduce the TQ relation, though in 
this case it can be done very easily (see, for example the recent ODE/IM review 
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article [3]). Correspondingly, x becomes subdominant at a: ^ oo on the positive 
real axis: xi^i^iJ) ^ ij{x,Ei,l). The motivation of dealing with instead of 
ip (jl.6p . is two-fold. Firstly, x can be obtained by applying the powerful and 
numerically efficient iterative method proposed by Cheng many years ago [9] in 
the context of Regge pole theory, and applied to spectral problems of this sort in 
|,10]. To this end we introduce the linear operator L, defined through its formal 
action 

LlxP] = -7 ; . (1.14) 

^ ^ {p+l){p-l-l) ^ ' 

So for any polynomial V{x) of x , 



^--](^ + -]L[V{x)\=V{x), (1.15) 
ax X \ax X 



and the basic differential equation (jl.Sp . with the boundary conditions (|1.9p at 
the origin, is equivalent to 

X(x, E, I) = x'+^ + L [{x^^' - E)xix, E, I)] . (1.16) 

Equation (jl.l6p is solvable by iteration and it allows the predictions of the 
ODE/IM correspondence to be checked with very high precision. 

The initial results of [H [21 [3] connected conformal field theories associated with 
the Lie algebra Ai to (second-order) ordinary diff'erential equations. The gener- 
alisation to ^n-i-models was established in [HI [12] but it was only recently [l3] 
that the ODE/IM correspondence was generalised to the remaining classical Lie 
algebras Bn, Cn and D„. Our attempts to derive generalised TQ relations from 
the proposed set of pseudo-differential equations were unsuccessful, but a series 
of well-motivated conjectures led us directly to the BAE, allowing us to estab- 
lish the relationship between BAE and pseudo-differential equation parameters. 
Moreover, while the numerics to calculate the analogs of the functions ip turned 
out to be very costly in CPU time, the generalisation of Cheng's method proved 
very efficient and allowed very high precision tests to be performed. This is our 
second main reason to deal with solutions defined through the behaviour about 
X = 0, rather than x = 00. 



2 Bethe ansatz for classical Lie algebras 

For any classical Lie algebra g, conformal field theory Bethe ansatz equations 
depending on a set of rank{Q) twist parameters 7={7a} can be written in a 
compact form as 

rank{s) ^(b) f rp{ci) n 

rr n^a.,.^B^S^i_2L = _-^^ z = 0,1,2,... (2.1) 
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where q'^^\e,^) = Q^'^\Qf^E,^), and the numbers i?!"^ are the (in general com- 
plex) zeros of the functions Q^"'\ In (j2.ip the indices a and b label the simple 
roots of the Lie algebra g, and 

Bab= , ^""'"'^ 12 , a,6 = l,2,...,ranM0) (2.2) 
I long roots I ^ 

where the a's are the simple roots of g. The constant J7 = exp {i-^-ji^ is a pure 
phase, /X is a positive real number and /i^ is the dual Coxeter number. 

It turns out that the Bethe ansatz roots generally split into multiplets (strings) 
with approximately equal modulus |-E'i"''|. The ground state of the model corre- 
sponds to a configuration of roots containing only multiplets with a common 
dimension da = K/Baa] the model-dependent integer K corresponds to the de- 
gree of fusion (see for example |14]). 

3 The pseudo-differential equations 

To describe the pseudo-differential equations corresponding to the An-i, Bn, Cn 
and Dn simple Lie algebras we first introduce some notation. We need an n^^- 
order differential operator [12] 



Dnig) = D{gn-i - (n-1)) D{gn-2 - (n-2)) . . . D{gi - 1) Dig^) , (3.1) 

depending on n parameters 

g={5n-i, • • • ,c/i,5o} , ^"={1^-1- go,n-l- gi,...,n-l- gn-i} ■ (3.3) 

Also, we introduce an inverse differential operator {d/dx)~^^ generally defined 
through its formal action 

-r\ = , (3.4) 

dx) s + l ^ ' 

and we replace the simple 'potential' P{E, x) = (x^*^ — E) of (II. Sp with 

Pk{E,x) = {x^"'^I'' -Ef . (3.5) 



Using the notation of Appendix B in [13] the proposed pseudo-differential 
equations are reported below. 

An-i models: 

The An-i ordinary differential equations are 

Dr,{g^)xi_^{x,E) = PK{x,E)xi_i{x,E) , (3.6) 
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with the constraint Y17=o 9i= ^^\ ™d the ordering gt < gj < n — 1, \/ i < j . 
We introduce the alternative set of parameters 7 = 7(g) = {7a (g)} 

^^ = Wm[Es^-^^) . (3.7) 

The solution xli-ii^^ E) is specified by its x ~ behaviour 

^ x"^^^^" + subdominant terms, {x 0"*") . (3.8) 

In general, this function grows exponentially as x tends to infinity on the positive 
real axis. In Appendix B of |13j . it was shown that the coefficient in front of 
the leading term, but for an irrelevant overall constant, is precisely the function 
Q^^\E,^) appearing in the Bethe Ansatz, that is 

Xn-i ~ Q 2 e M+i + subdominant terms, {x ^ 00) . (3.9) 

Therefore, the set of Bethe ansatz roots 

{i^f)}^Q«(i?«,7) = (3.10) 
coincide with the discrete set of E values in (13.61) such that 



xl-i^ oix'^^-'^'-e'^w \ , (x^oo). (3.11) 
This condition is equivalent to the requirement of absolute integrability of 

/^(n-i)-e-M+rj^t_^(^,i^) (3.12) 

on the interval [0, 00). It is important to stress that the boundary problem defined 
above for the function Xn-i (|3.8p is in general different from the one discussed 
in Sections 3 and 4 in [13j involving ^{x^E). The latter function is instead a 
solution to the adjoint equation of ()3.6p and characterised by recessive behaviour 
at infinity. Surprisingly, the two problems are spectrally equivalent and lead to 
identical sets of Bethe ansatz roots. 

Dn models: 

The Dn pseudo-differential equations are 

^n(g^) ^n(g)X2n-l(x, E) = ^Pk{x,E) {J-^ ^/Pk{x,E) X2n-l{x, E) . 

(3.13) 
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Figure 1: Lowest three functions ^{x,E) for a D4 pseudo-differential equation. 



Fixing the ordering gi < gj < h'^ /2, the g <-> 7 relationship is 



la 



2K 

Wm 



Y.di-^h''] , (a = l,...,n-2) 



(3.14) 



7n-l 



K 



^9i- ^h"^ , In 



Z^9i- 9n-i 



,4=0 / \i=0 

The solution is specified by requiring 

X2n-i ~ x'^^"^" + subdominant terms, {x — > 0"*") 



(3.15) 



(3.16) 



X2n-i ~ Q (-E',7(g)) ~e + subdominant terms, {x —>■ 00). (3.17) 



Figure [H illustrates ^{x,E) = x'^^ 2 e J^^+i X2n~iix,E) for the first three eigen- 
values of the 1)4 pseudo-differential equation defined by K=1,M = 1/3 and 
g=(2.95,2.3,l.l,0.2) . 

Bn models: 

The Bn ODEs are 

Dnig^)Dnig)xL-ii^,E) = VPk{x,E) (J-^ ^PKix,E)xL-iix,E) . (3.18) 



With the ordering g^ < gj < /i^/2, the g ^ 7 relation is 

2K ("-^ ^ 

la 



.4 = / 



(3.19) 



The asymptotic behaviours about x = and x = 00 are respectively 



X2„_i ^ x^ ^° -|- subdominant terms, (x — > 0"^) , 



(3.20) 
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and 



xln-i ~ Q^'^H^nis)) X ^" 2 e M+i + subdominant terms, ^ oo) . (3.21) 
C„ models: 

The pseudo-differential equations associated to the C„ systems are 

Dn{^^) D„(g)x2n+l(x,^;) =Pi^(x,£;)(^^^ Pk{x , E) X2n+l{x , E) 

(3.22) 

with the ordering gi < gj < n. The relation between the ^'s and the twist 
parameters in the BAE is 

> = (g 9^ - ' 7n = ^ (^g 5. - (3.23) 

and 

xL+i ~ x^""^'' + subdominant terms, (x ^ 0+) , (3.24) 

xL+1 ~ Q'^^H-E', 7)x"''*^e^^^ + subdominant terms, (x ^ oo) . (3.25) 



Using a generalisation of Cheng's algorithm, the zeros of Q^^\E,'y) can be 
found numerically and shown to match the appropriate Bethe ansatz roots |13j . 

In general, the 'spectrum' of a pseudo-differential equation may be either 
real or complex. In the An^i, Bn, Dn models with i^=]0, the special choice 
Qi = i leads to pseudo-differential equations with real spectra, a property which 
is expected to hold for a range of the parameters g (see, for example, [E]). The 
K>1 generalisation of the potential (13. 5p . proposed initially by Lukyanov for the 
Ai models [15] but expected to work for all models, introduces a new feature. 
The eigenvalues corresponding to a K=2, 3 and iC=4 case of the SU{2) ODE 
are illustrated in figure [2j The interesting feature appears if we instead plot the 
logarithm of the eigenvalues as in figure [3l We see that the logarithm of the 
eigenvalues form 'strings', a well-known feature of integrable models. The string 
solutions approximately lie along lines in the complex plane, the deviations away 
from which can be calculated [13j using either WKB techniques, or by studying 
the asymptotics of the Bethe ansatz equations directly. 



* The Cn spectrum is complex for any integer K>1. 
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Figure 2: Complex ii^-plane: the eigenvalues for the SU (2) model with M = 3, 
5o = for = 2, 3 and 4 respectively. 
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Figure 3: Complex (In £')-plane: two, three- and four-strings. 



To end this section, we would like to comment briefly on the motivation behind 
the conjectured pseudo-difi'erential equations of Bn, Cn and Dn type. Modulo the 
generalisation to K>1, the An-i type ODEs were derived in [12]. We began with 
the D3 case since it coincides up to relabeling with ^43, implying that the -D3 
function Q^^\E,^) coincides with the function Q^'^\E,'y). Fortunately, the 
latter is known [12] to encode the spectrum of a differential equation satisfied by 
the Wronskian of two solutions of the Q^^-'-related ODE. The generalisation to 
Dn models with larger n was then clear. Further supporting evidence came from 
a relationship between certain Z?„ lattice models and the sine-Gordon model, 
which appears as an SU{2) problem. This relationship also extends to a set of 
Bn models, and leads naturally to the full Bn proposal. Finally, the C„ proposal 
arose from the Bn cases via a consideration of negative-dimension W-algebra 
dualities [TB]. Numerical and analytical tests provided further evidence for the 
connection between these spectral problems and the Bethe ansatz equations for 
the classical Lie algebras. 

4 Conclusions 

The link between integrable models and the theory of ordinary differential equa- 
tions is an exciting mathematical fact that has the potential to influence the 
future development of integrable models and conformal fleld theory, as well as 
some branches of classical and modern mathematics. Perhaps the most surpris- 
ing aspect of the functions Q and T, only briefly discussed in this short note, is 
their variety of possible interpretations: transfer matrix eigenvalues of integrable 



lattice models in their CFT limit OE], spectral determinants of Hermitian and 
PT-symmetric \n\ [TH] spectral problems (see for example pLOj), g- functions of 
CFTs perturbed by relevant boundary operators [5l [19], and particular expec- 
tation values in the quantum problem of a Brownian particle |20) . Further, the 
(adjoint of the) operators (j3.6p . (|3.13p . (|3.18p and (j3.22p resemble in form the 
Miura-transformed Lax operators introduced by Drinfel'd and Sokolov in the 
context of generalised KdV equations, studied more recently in relation to the 
geometric Langlands correspondence |21l [22] . Clarifying this connection is an in- 
teresting open task. Here we finally observe that the proposed equations respect 
the well-known Lie algebras relations D2 ^ Ai(BAi, D^, Bi ~ ^i, i?2 ~ C'2- 
Also, at special values of the parameters the Cn equations are formally related 
to the Dn ones by the analytic continuation n —n, matching an interesting 
W-algebra duality discussed by Hornfeck in (16] : 

{D-n)K X {D-n)L ^ {Cn)-K/2 X (C'n)-L/2 ^g) 
{D_n)K+L iCn)-K/2-L/2 

The relationship between our equations and coset conformal field theories is an- 
other aspect worth investigation. We shall return to this point in a forthcoming 
publication. 
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